The J- and H-bands of dye aggregate spectra: Analysis of the coherent exciton 

scattering (CES) approximation 
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The validity of the CES approximation is investigated by comparison with direct diagonalisation 
of a model vibronic Hamiltonian of TV identical monomers interacting electronically. Even for quite 
short aggregates (N > 6) the CES approximation is shown to give results in agreement with direct 
diagonalisation, for all coupling strengths, except that of intermediate positive coupling (the H-band 
region). However, previously excellent agreement of CES calculations and measured spectra in the 
H-band region was obtained [A. Eisfeld, J. S. Briggs, Chem. Phys. 324 (2006) 376} 1 . This is shown 
to arise from use of the measured monomer spectrum which includes implicitly dissipative effects 
not present in the model calculation. 
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I. INTRODUCTION 

In a scries of papers^— the "coherent exciton scatter- 
ing" (CES) approximation was applied with some success 
to calculate the absorption and circular dichroic spec- 
tra of a variety of aggregates of dye molecules. Accord- 
ing to their structure, such aggregates exhibit continu- 
ous J- or H-band absorption (see e.g. Refs££) and the 
CES approximation, working directly with the continu- 
ous vibronic absorption profile, is able to reproduce both 
types of bandshape in detail. The CES approximation 
includes intra-molecular vibrations explicitly and cou- 
pling to other modes only implicitly by working with the 
continuously-broadened experimental spectra. Over the 
years there have been many other approaches to the in- 
clusion of vibronic effects on aggregate spectra, see for 
example Refsji— . 

The CES approximation is a mean-field type of ap- 
proximation in which the exact monomer Green func- 
tion, involving both electronic and vibrational degrees 
of freedom, is replaced by its average in the vibrational 
ground state. This allows the aggregate bandshape to 
be expressed in terms of the monomer absorption band- 
shape; the only fit parameter being the inter-monomer 
electronic interaction V. If G(E) = (E - H M -V + iS)- 1 
is the aggregate Green operator at energy E and g(E) = 
(E - H M + iS)- 1 that of the non-interacting monomers, 
one has the identity 



G 



gVG 



(1) 



where Hm is the sum of monomer Hamiltonians and V is 
the electronic interaction operator between monomers. 
The CES approximation corresponds to replacing the 
exact monomer Green operator g by (g), where (...) 
denotes the aggregate vibrational ground state average. 
Then 



becomes 



(G) = (g) + (gVG) 



(G) = (g) + (g)V(G) 



(2) 



(3) 



when V is assumed to be independent of vibrations. 
Eq. ([3]) has symbolically the solution 



(G) 



i - (g)v 



<<?>• 



(4) 



However, the physical content of the CES approximation 
is best seen by iterating Eq. ([3]) in a Born series, i.e. 



(G) = (g) + (g)V(g) + (g)V(g)V(g) 



(5) 



If the matrix element Gr 



(-K n \G\-K m ), where \n n ) is a 



state in which monomer n is excited electronically with 
all other monomers in their ground electronic state, is 
taken, then Eq. ([5]) becomes 



(G r , 



:(gn)Snm + (gn)V nm (g m ) 



(gn)V nn > (gn')Vn' 



(6) 



This shows that whenever electronic excitation is handed 
on from one monomer to another, the monomer which de- 
excites goes back into its ground vibrational state. This 
is the essential ingredient in the CES approximation. 

The CES approximation has been used successfully to 
describe the detailed lineshape of various measured ag- 
gregate spectral— in the intermediate and strong cou- 
pling regime. Following Simpson and Peterson^ the 
strong and weak coupling regimes are defined by the con- 
ditions V > A or V < A respectively, where A is the 
width of the monomer absorption spectrum. One speaks 
about intermediate coupling when V ~ A. 

In some ways, the CES approximation performs better 
than it ought to, in that in its original formulatio n 17 i 18 
it involves only the (zero temperature) aggregate ground 
state, represented as a simple product of monomer vi- 
brational ground states. However in the cases to which 
it has been applied, measurements were made at room 
temperature in a variety of solvents which also inter- 
act with the dye molecules. The origin of the success 
of CES in this more general context may lie in the recent 
demonstration^ that its essential structure (Eq. (Q}) is 
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retained when the average (. . .) is extended to include 
finite temperature and averages over site solvent interac- 
tions. This question of temperature and solvent effects 
on J-aggregate spectra has been studied extensively both 
experimentally^ and theoretically^. In the theoretical 
work a purely electronic Hamiltonian was used for the 
aggregate but coupling to phonons of the surroundings 
were considered to show that temperature effects can 
be included successfully in this model. Here, since the 
CES method includes such effects implicitly, coupling to 
phonons of the surrounding will not be considered explic- 
itly. Rather, we will include the intra-monomer vibra- 
tions and compare the results of the CES approximation 
with those of a direct diagonalisation of the aggregate 
vibronic Hamiltonian. 

Although the original derivatio n 17 ' 18 suggests that the 
CES-approximation is applicable to strong coupling, in 
Ref^ it was claimed that it should only be valid for weak 
coupling. Our subsequent comparison with experiment 1 
indicates validity for both strong and intermediate cou- 
pling. To further throw light on these questions, in this 
work the validity of the CES approximation will be exam- 
ined in detail in its original form, in which an average is 
taken over the vibrational ground state only. In sectionllll 
the method of calculation is described. In section UTT1 the 
results are presented and in section IIVI our conclusions 
as to the limits of validity of the CES method are given. 



n is in the excited electronic state and all others in the 
ground state, as 



Kn) = lxi)IX2> • ■ ■ \<t> n ) ■ ■ ■ \xn)- 



(8) 



Here we generalise this notation to define an aggregate 
vibronic state in which monomer n is excited electroni- 
cally and in vibrational state |/3„) and all other monomers 
in the ground electronic state with vibrational quantum 
number a, for monomer i. This generalised basis is writ- 
ten 



«1 



\Pn) 



\OtN) 



(9) 



i.e. u n denotes the sequence of vibrational quantum 
numbers (a\ . . . /3 n . . . ajv)- Note that this simple prod- 
uct basis is a BO basis for the aggregate as a whole. 
However, the inter-monomer interaction mixes these ba- 
sis states to give no BO separation of the total aggregate 
eigenfunctions. In this chosen basis the matrix elements 
of the aggregate Hamiltonian are 



N 



+ V nm ({3 n \a n )(a m \fl m ) S ai 



i=l 



(10) 



II. THEORETICAL METHOD FOR 
CALCULATION OF AGGREGATE SPECTRA 



with 



(11) 



Direct diagonalisation of the vibronic 
Hamiltonian 



We will consider an aggregate consisting of N identical 
monomers having only one excited electronic state upon 
which vibrational manifolds are built. For an arbitrary 
monomer n, the ground state is the product |Xn)l a n) 
where \xn) is the electronic ground state and \a n ) de- 
notes the vibrational state with quantum number a n in 
this electronic state (we assume Born-Oppenheimer (BO) 
separability within a monomer). In a corresponding no- 
tation the excited vibronic state of monomer n is desig- 
nated by the product \4> n )\/3 n ). Throughout this work x 
and a will be used for the electronic ground state and its 
vibrational state and <p an d P for the excited electronic 
state and its vibrational state. 

The aggregate Hamiltonian is the sum of the monomer 
Hamiltonians H n and the inter-monomer interaction op- 
erator V 



H A = H M + V = J2 H n + V. 



(7) 



In the following it is assumed that the electron-electron 
interaction V is independent of nuclear co-ordinates. In 
previous worker— i 17 i 18 (and as used in Eq. ©) we defined 
one-exciton aggregate electronic states, where monomer 



Here e' 3 ™ is the energy of monomer n in the vibrational 
state \j3 n ) of the excited electronic state and e ai that of 
monomer i in the vibrational state |a,) of the electronic 
ground state. The quantum numbers a n and (3 m are 
elements of the sequence u m . The Franck-Condon (FC) 
factor (a n \(3 n ) denotes the overlap of vibrational state /3 n 
of the excited electronic state with vibrational state a n 
of the ground electronic state. 

The eigenstates and eigenenergies of the aggregate are 
defined by 



(12) 



We designate by \n, 0) the absolute vibronic aggregate 
ground state i.e. 



ko> = nix*>k- = >- 



(13) 



Then the aggregate absorption cross-section energy de- 
pendence is given by the absorption strength 



A(Ei)= (7T,0|£-£> A |tfi) 



(14) 



where £ is the light polarisation vector and the aggregate 
dipole operator Da is the sum of monomer operators 
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D n . To evaluate Eq. (|T4|) further an identity is inserted 
between Da and \^i) so that 



^2 S ■ (n 7 0\^2D n \Tr m ,v> m )(TT m ,v m \^i 



The dipole matrix element appearing in Eq. (|15|) can be 
simplified to give 



{■K,Q\D n \lT m ,U m ) 

= {a n =Q\{Xn\Dn\<t>n)\l3n) <5nm<^ m , (0. . ./3„ . . .0) • 



(16) 



Consistent with our assumption that V is independent of 
vibrations, we make the Condon approximation that D n 
is independent of vibrations so that the dipole transition 
matrix element of monomer n is 

{a n \(Xn\D n \<i> n )\/3 n ) = d n {a n \/3 n ). (17) 

The absorption strength Eq. (fT4")l then reduces to 

= 0|/3„)<7r n ,(0... ft... 0)|*j 

Tl,Pn 

The aggregate (stick) spectrum is then given by 

A(E)=Y / AEi)S(E-E l ). (19) 
i 

The dimension of the basis Eq. (j9]) and therefore the 
number of terms in the summation Eq. (I19|) depends on 
the number of vibrational states taken into account. If 
all monomers have n g vibrational states of the ground 
electronic state and n e vibrational states of the excited 
electronic state, the dimension is given by 



dim = N ■ n e 



,iV-l 



(20) 



Therefore the number of absorption sticks is already large 
(about 10000) for an aggregate of ten monomers with 
two vibrational states in each electronic state. Note that 
the dimension depends crucially on the number n g of 
vibrational states of the electronic ground state. 

In the special case that we assume that all N monomers 
are identical and arranged in a circle, the dimension of 
the matrix to be diagonalised can be reduced by a fac- 
tor N using a transition to a delocalised excitonic ba- 
sis rather than the localised basis Eq. ©. For sim- 
plicity we specialise V to nearest-neighbour coupling. 
How the delocalised states are constructed is illustrated 
first by considering the simple case of a cyclic trimer 
(iV = 3). The excitonic basis consists of states \k,v>i}, 
with k = 27r/3, 47r/3, 27r, which are simply a superposi- 
tion 



=-7= [e lfe |7ri,(/3i,"2,a 3 )) 
V iV 

+ e 2ifc |7r 2 ,(a 3 ,/3i,a 2 )) 
+e 3ik \n 3 ,(a 2 ,a 3 , ft))] 



(21) 
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FIG. 1: Sketch of the action of the shift operator S. As ex- 
ample 52 is considered. The excited electronic state is high- 
lighted in gray. 



of the three possible shifts of the |7Ti, Vy) state along the 
cyclic aggregate. In general, for a cyclic iV-mer the exci- 
tonic basis \k,u{) is 



N 
n=l 



|fc )I / 1 > = ^=>> 8fc "S n _ 1 | 7 r 1 , I / 1/ 
V jv 



(22) 



with k — 2irj/N, j — 1, . . . , N and where S m is a shift 
operator shifting both the electronic and vibrational ex- 
citation on the aggregate by m monomers to a higher 
monomer index i ' = i + m, i.e. 

)>■ (23) 

The action of S m is also illustrated in Fig. Q] 

Denoting by W = V n , n+ i the nearest-neighbour inter- 
action between the monomers and defining 



C 



W for N — 2 (dimer) 
2W for N > 2 



(24) 



the matrix elements of H A in the exciton basis Eq. 
are 



(^ialI-HaU.^i) = 



,C 



N-l 



c lk -{f3 1 \a 2 ){a N \Pi) ]J 



.c 



i=2 
N-l 



e lk -{Pi\a N ){a 2 \h) ]J S ai+li 



i=2 



8 k j- 
(25) 



Note that Ui = (ft, a 2 , ■ ■ ■ , a/v). From Eq. (|25)l one 
sees that Ha is diagonal in the exciton k basis i.e. the 
eigenvalue problem 



k \ 



Ekxl^kx) 



(26) 



can be solved separately for each k = 2w/N, . . . , 2tt with 
A = 1, . . . , n e ■ n^ -1 (n g and n e are the numbers of vi- 
brational states taken into account in the ground and 
excited electronic state of each monomer, see Eq. ([20|) ). 
The absorption strength (see Eq. (US])) from the aggre- 
gate ground state \ir, 0) to the eigenstate \^k\) is given 
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by 



1 i N 



71=1 

n.-l 



5^(ai = 0|A)<*,( j 8i > 0...0)|*j kA ; 



3i=0 



(27) 



The first factor on the r.h.s. is the effective electronic 
dipole moment of the aggregate in exciton state k and 
depends upon the aggregate geometry. In the simplest 
case where all monomer dipoles are parallel, this factor 



reduces to N 



e-di 



Sk,2TT- In this case, to calculate the 
absorption spectrum, only the block k — 2n of the Hamil- 
tonian matrix (j2"5]l needs to be diagonalised and the stick 
spectrum is given by (see Eq. (1191) ) 



(28) 



A(E) = A(E 27t , x )5(E - E 2 ^ x ) 



with 



A(E 2ViX ) = N E-d^ («i=0|/?i) (k, (A, ... 0)|* 



2w,X/ 



(29) 



Calculations of the aggregate spectrum from Eq. (j29j) us- 
ing the eigenfunctions and eigenenergies from direct di- 
agonalisation of the aggregate Hamiltonian will be pre- 
sented and referred to as DD calculations. 

In the following discussion the absorption strength is 
given in arbitrary units. All energies, including a and A 
are given in units of Tuo. 



B. The monomer spectrum 

All that is required to specify the vibrational shape of 
the monomer spectrum and carry through the calcula- 
tion of the aggregate spectrum are the monomer Franck- 
Condon factors between vibrational states of upper and 
lower electronic states (for the definition see the discus- 
sion after Eq. (fTTj) ). Hence the calculations can be made 
for arbitrary ground and excited state BO potential sur- 
faces. However, to conform with a large body of work, 
particularly using the second quantisation formalism, we 
will consider here that upper and lower potentials are 
both harmonic of the same curvature, with a shift Q of 
the minimum of the upper surface w.r.t. the minimum 
of the lower. Here Q is the vibrational coordinate giv- 
ing rise to a single dominant vibrational progression in 
the monomer spectrum. Despite its simplicity this model 
gives a surprisingly good fit to the spectra of a large va- 
riety of organic molecules. It has the advantage that 
the FC factors can be calculated analytically^ and the 
monomer spectrum obtained in a simple analytic form. 
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FIG. 2: Monomer spectra of different width (standard de- 
viation) A convoluted with narrow Gaussians (of standard 
deviation a = 0.1). All energies are in units of foj. 



Denoting the harmonic vibration frequency by w and 
defining the Huang-Rhys factor— by 



X 



— Q 2 



(30) 



the absorption strength of the monomer from the ground 
state (where a = 0), which is of relevance to the CES 
approximation, is 



M{E) 



n c -l 

E 

/3=0 



-x 



0\ 



-S(E- pfwj). 



(31) 



The width (standard deviation) of the monomer spec- 
trum in this case of a Poissonian distribution is given 
by 



A 



xhu. 



(32) 



Examples of monomer spectra for various values of the 
parameter A are shown in Fig. [2] As is well-known, for 
small X , i.e. small A, the Poisson distribution gives only 
one dominant peak (Fig.^a)) but becomes progressively 
more Gaussian as X increases and the progression broad- 
ens. 



C. The CES approximation 

It is clear from the discussion of the introduction that 
the CES approximation is equivalent to performing the 
diagonalisation of section III Al in a basis in which in the 
ground electronic state only the ground vibrational state 
is included. Then transfer of excitation to or from a 
monomer must proceed through this state as can be seen 
from the expansion ^ of the CES approximation. How- 
ever the effect of the CES approximation on the absorp- 
tion lineshape is better illustrated by the form (|4|. The 
imaginary part of the monomer function (g(E)) is pro- 
portional to the monomer absorption strength Eq. pip . 
Once Im(g) = g 1 is determined, the real part can be 
calculated from the Kramers-Kronig dispersion relation. 
Then the aggregate absorption strength can be calculated 
from Im(G), where (G) is given by expression 
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FIG. 3: Monomer function 1/<?h as a function of energy 
E. The width of the monomer spectrum is A = 0.8. The 
absorption strength at each pole is indicated by the size of the 
circle. The dashed horizontal lines are for energies C = —1.5 
and C = 1.2. 
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FIG. 4: Illustration of the convergence criterion for the ag- 



gregate absorption. 
The overlap is 68%. 
The overlap is 81%. 
The overlap is 90%. 
The overlap is 95% 
A = 0.8, n e = 9. 



a) n g = 1 (solid) and n g = 2 (dashed). 

b) n g — 2 (solid) and n g = 3 (dashed). 

c) n 9 = 3 (solid) and n g = 4 (dashed), 
b) n s = 4 (solid) and n g = 5 (dashed). 
Used parameters are JV = 4, C — 1.2, 



III. NUMERICAL RESULTS 



For a stick monomer spectrum, the Re(g) = is read- 
ily calculated and the function 1 / g R (E) is shown in Fig. [3] 
(see also Ref»i). This function diverges at the discrete 
energy points where the monomer absorbs (given by the 
energy values /3Hui). Outside of these points the function 
(g) is wholly real so that the aggregate function (G) has 
poles (and hence there is absorption) where the coupling 
strength C — l/g R (for the definition of C see Eq. pi])). 
These poles are given by the intersection of horizontal 
lines (representing different coupling strength C) with 
the function l/g a (E). From Fig. [3] one sees that the 
location of aggregate absorption depends sensitively on 
the sign and magnitude of C . As discussed in detail in 
Refi this explains the qualitative dissimilarity of J and 
H band spectra. At the aggregate absorption energies 
(poles) Ei, g, — so the expression (|4]) is not defined. 
Nevertheless^ one can calculate the absorption function 
of the aggregate as 



This shows that the absorption strength at a particular 
pole is inversely proportional to the slope of g" 1 at that 
point. Hence only the poles where g" 1 has gradient in 
the order of unity absorb appreciably and these points 
are denoted by circles in Fig. [31 One also notes that in 
the CES approximation the number of poles is limited by 
the number n e of vibrational states included on the up- 
per potential irrespective of the number N of monomers 
considered. 



All calculations, both direct diagonalisation (DD) and 
CES calculations yield stick spectra. The CES results in 
fact were obtained from the DD calculation by restricting 
to the ground vibrational state only (i.e. n g = 1), but it 
was checked that the results agreed with those obtained 
from Eq. ((33)) . To facilitate comparison of CES and DD 
results the calculated stick spectra have been convoluted 
with a narrow Gaussian of standard deviation a = 0.1 tbu. 
As can be seen in Fig. [2] this width is small enough to 
clearly resolve the individual peaks of the monomer ab- 
sorption spectrum. 

The dimension of the Hamiltonian matrix (see 
Eq. (|25p ) is n e ■ n^ -1 and computer time for diagonal- 
isation therefore escalates, particularly as n g increases. 
However all DD spectra have been checked for conver- 
gence and the cases presented have a minimum 95 % 
overlap with the spectra with n g and/or n e increased 
by unity. For all spectra shown, we have used n e = 9. 
In Fig. [5] our convergence criterion is illustrated for an 
typical example. 

There are three parameters deciding the shape of the 
aggregate spectrum, namely, the coupling strength C, 
the monomer spectral width A and the size N of the 
aggregate. In the comparisons presented, representative 
points in this 3-dimensional space are chosen. To some 
extent parameters C and A can be combined into the 
Simpson-Peterson (SP) parameter |C|/A, where strong 
coupling has |C|/A > 1 and weak coupling |C|/A < 1. 
Although this is a very rough measure, it will emerge 
that, for N ^> 2, the CES approximation is good for 
both weak and strong coupling and only may have prob- 
lems in the intermediate positive regime. Similarly, it 
has been shown that the CES approximation is not good 
for dimers (N = 2), since electronic excitation is trapped 
indefinitely on just two monomers (see Refj2i). This con- 
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elusion is reinforced by the results presented here, except 
in the case of weak coupling. 

In considering these comparisons, two sum rules sat- 
isfied by both DD and CES spectr a 17 i 25 should be kept 
in mind. One is that the total absorption strength of 
any spectrum is equal to that of the monomer (first sum 
rule SRI) and secondly that the mean energy of the spec- 
trum (first moment) is shifted from that of the monomer 
by exactly the coupling energy C (second sum rule SR2). 

Since the CES result is independent of N, it is logi- 
cal to present comparisons of DD spectra with CES for 
fixed values of C and A and see at which IV the CES 
approximation becomes adequate. 

Although we have performed calculations for monomer 
spectra with X, defined in Eq. (|30|) . ranging from X = 
to X = 2, in this work we will focus on the case X = 0.64, 
i.e. A = 0.8hu. This value of X is representative of 
those found for the dominant vibrational progression of 
many organic dyes. 

1) Strong negative coupling 

This case is the classic case of the formation of a J 
band. Representative spectra are shown in Fig. [5] (for 
SP fa 6). The values of n g indicated on the figures are 
those used in the calculation of the DD spectra. The 
results presented are for C — — 5 but we have checked 
that there is no change in the shape of the spectra 
for C < — 5. The monomer spectrum for this case is 
given in Fig. [He). As N increases, essentially only one 
strong peak on the low-energy side of the monomer band 
appears, whose mean energy is shifted by C from the 
monomer mean energy. In Fig. 03(d) (N — 8) one sees 
that the CES approximation is in good agreement with 
the DD result. As a result of the sum rules, the small 
peak near E — —3.5 in the DD spectrum causes its J 
band to be slightly less intense and shifted to slightly 
larger negative energy than the CES J band peak. On 
the dispersion curve of Fig. [3] one sees the origin of the 
J band as due to the single pole at negative energy split 
off from the monomer band. This pole carries all the 
oscillator strength. For the other poles the g~ x curve is 
almost vertical, giving vanishing absorption according 
to Eq. ([33|. 

2) Intermediate negative coupling 

Here we choose C = —1.5 so that SP ~ 2, all other 
parameters being as in Fig. [5j From Fig. [5] one notes 
that now the CES approximation has small subsidiary 
peaks. As can be seen from Fig. [3] these are due to 
additional poles in the monomer region acquiring some 
oscillator strength. Nevertheless the CES spectrum is 
dominated by the isolated-pole J-band peak. Although 
the peaks in the region of monomer absorption are 
different in the DD calculation, these contributions 
diminish with increasing N such that for N = 8 there is 
again very good agreement between CES and DD results. 
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(a) N = 2 

n g = 8 


Ill 


(b) N = 4 

n g = 8 




(c) N = 6 

n g = 4 




(d) N = 8 

n g = 2 



-5 -4 -3 -2 

Energy [huj] 

FIG. 5: Strong negative coupling. Full curves CES spectra, 
dashed curves DD spectra, for C = —5, A = 0.8. The values 
of TV and n g are indicated on the figures. 



3) Weak negative and weak positive coupling 

Here C = -0.2, SP = 0.25 (Fig. EJ) and C = +0.3, 
SP pa 0.38 (Fig. all other parameters being taken as 
in Fig. [5j 

In the case of weak coupling the DD spectrum does 
not change with IV for the case A = 0.8 and so only 
the dimer case N — 2 is shown in Fig. [7] and Fig. [5] 
Here there is excellent agreement, both in position and 
magnitude between DD and CES results. This is perhaps 
remarkable when one remembers that the DD results rep- 
resent a convolution over up to about 10000 individual 
eigenvalues, whereas the CES spectrum arises from only 
the four main contributing poles for C = —0.2 and for 
C = +0.3, respectively. 

To summarise, we have shown that, for aggregates of 
length N > 6, the CES results, i.e. including only the 
ground vibrational state in the electronic ground state, 
agree well with DD calculations where several vibrational 
states are included, for all negative coupling strength and 
for weak positive coupling. For weak negative or positive 
coupling the CES method gives good results even for the 
dimer N = 2. 

The case of strong positive coupling is not shown, since 
if there were a strong shift to higher energies the absorp- 
tion would overlap higher electronic bands and our model 
would break down. 

There remains the most interesting case of interme- 



7 



3 
-ft 

& 6 

a 

o 

+3 3 
ft 

o 

,0 



(a) TV = 2 



(b) JV = 4 



(c) TV = 6 



(d) TV = 8 



ft 



.ft 




-10 12 

Energy [ftoj] 



FIG. 6: Intermediate negative coupling. Full curves CES 
spectra, dashed curves DD spectra, for C = —1.5, A = 0.8. 
The values of TV and n g are indicated on the figures. 
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FIG. 8: Weak positive coupling. Full curve CES spectrum, 
dashed curve DD spectrum, for C = +0.3 and A = 0.8 with 
TV = 2 and n„ = 8. 
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FIG. 7: Weak negative coupling. Full curve CES spectrum, 
dashed curve DD spectrum, for C = —0.2, A = 0.8, with 
TV = 2 and n„ = 8. 



FIG. 9: Intermediate positive coupling. Full curves CES spec- 
tra, dashed curves DD spectra, for C = +1.2, A = 0.8, convo- 
luted with a narrow Gaussian of standard deviation o — 0.1. 



diate positive coupling where the aggregate absorption 
occurs in the region of monomer absorption, typical of 
broad H band formation. This case will be considered in 
some detail. 

4) Intermediate positive coupling 

Note first that from Fig. [31 one sees that for interme- 
diate positive coupling, the aggregate poles in the CES 
approximation lie in the region of monomer absorption 
and their strength and location then is very sensitive to 
the coupling strength C. In the following we consider 
the case C = 1.2 (i.e. a SP parameter of 1.5). In Fig. [3J 
the CES result for this case is illustrated (intersection of 
upper dashed horizontal line) . In Fig. [S] a comparison of 
CES and DD spectra is shown. 

One notes no good agreement of CES and DD results 



in this case, except to say that the dominant absorption 
is in the same spectral region for both methods. Indeed 
in this case the DD results themselves vary significantly 
with the TV value and no convergence with TV is seen, 
in contrast to the case of intermediate negative coupling 
shown in Fig. |B1 

One must remember that the curves of Fig. [S] are cal- 
culated, in all cases, first as stick spectra and then convo- 
luted with a Gaussian whose width is chosen arbitrarily 
as (7 = 0.1. This width is considerably less than that 
typical of room-temperature solvent spectra. Hence in 
Fig. [10] the same data as in Fig. [9] is shown convoluted 
with a Gaussian of a = 0.3 width. Then one sees that, 
for TV > 4, there is a much better, although still not 
exact agreement between CES and DD results. Hence, 
even when averaged over a broad energy region, the CES 
method would appear not to perform well for the H- 
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FIG. 10: Same as Fig.[9]but convoluted with a broad Gaussian 
with a — 0.3. 



band case of intermediate positive coupling. However, 
in Ref.— we showed that the CES method gives excellent 
agreement, even in some fine details, with measured ri- 
band spectra. Important is that in Ref. 1 the measured 
monomer spectrum, rather than a calculated Poissonian 
stick spectrum, was used as input to the CES method. 
As we show now, it appears that this is crucial in order to 
obtain agreement with measured aggregate spectra and 
procedures based on fits to stick spectra in general appear 
less reliable. To illustrate this, specifically we adopt the 
following strategy of comparing three different methods 
of calculation of the aggregate spectrum. 

a) The measured monomer spectrum of Pinacyanol 
from Ref<2&, which was used in the CES calculations of 
Refi and is shown in Fig. UTTa) (dotted line), is fitted 
to a stick Poissonian distribution with each peak broad- 
ened by the same Gaussian. This fitted spectrum is also 
shown in Fig. [TTta) (solid line). Taking X and ui from 
the fitted monomer stick spectrum, the stick spectrum of 
the aggregate is calculated using the DD method. Then 
each peak is convoluted with the same Gaussian as used 
in the monomer spectrum. The coupling strength C has 
been optimised to obtain a best fit and the overall calcu- 
lated spectrum shifted by exactly the shift S = 480cm -1 
used in Refi for the dye concentration considered. The 
DD result for N = 5 and n g — 5 is shown in Fig. Efb) 
and is in reasonable agreement with experiment, except 
for a large discrepancy in the peak around 16, 500cm -1 . 

b) The fitted continuous monomer spectrum of 
Fig. IPTT a) is used as input in the CES calculation of the 
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FIG. 11: Dotted: measured Pinacyanol monomer (Fig. (a)) 
and aggregate spectrum (Fig. (b)-(d)) from Ref.— . a) Solid 
line: fit to a Poissonian convoluted with a Gaussian of width 
a — 450 cm -1 , b) Dashed: DD calculation with the peaks 
convoluted with the Gaussian of Fig. (a). In the calculation 
N = 5, n g — 5 was used, c) Solid line: CES spectrum ob- 
tained directly from the fitted monomer spectrum of Fig. (a), 
d) Solid line: CES spectrum obtained directly from the mea- 
sured monomer spectrum of Fig. (a). 



aggregate spectrum, again with C optimised and with 
S = 480 cm -1 . Now, Fig. ITTT c). there is reasonable 
agreement in the peak position and respective heights 
but the peaks appear too narrow. This probably arises 
from the too-low absorption of the monomer fit in the 
region around 19,000 cm -1 (see Fig. [TTta)). leading to 
a too-low absorption in the same region of the aggregate 
spectrum (as explained already, the local aggregate ab- 
sorption depends on the monomer absorption in the same 
locality). 

c) From Ref.—, the measured monomer spectrum is used 
as input to the CES calculation of the aggregate spec- 
trum. The result is shown in Fig.QjJd), where one notes 
excellent agreement with the measured aggregate spec- 
trum. 



IV. CONCLUSION 

In its original for m 17 ' 18 the CES approximation is 
equivalent to the assumption that each monomer pos- 
sesses only a ground vibrational state in the electronic 
ground state. Hence, in the aggregate only this state 
is involved when electronic excitation is transferred be- 
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tween monomers. Here the validity of this approximation 
has been tested by comparing CES spectra with those 
obtained from direct diagonalisation (DD) of the aggre- 
gate vibronic Hamiltonian with a sufficient number of vi- 
brational states in the electronic ground state to ensure 
convergence. For negative values of the inter-monomer 
coupling energy, leading to a shift of the aggregate ab- 
sorption to lower energies, good agreement of the CES 
spectra with DD spectra has been obtained for all nega- 
tive coupling strengths, so long as the aggregate size N 
exceeds 4 or 5 monomers. This success is probably due 
to the fact that for negative coupling, even intermedi- 
ate values, a single dominant J-band peak splits off from 
the monomer region and carries most of the oscillator 
strength. As shown in Fig. [3] this mechanism of J-band 
formation is well-described by the CES approximation. 

In the weak coupling regime, both negative and pos- 
itive, the CES result agrees excellently with DD results 
since there is no strong mixing of monomer vibronic tran- 
sitions. In fact for weak coupling CES and DD results 
agree even for the dimer N = 2. The success of CES in 
weak coupling can be understood as follows. For weak 
coupling one can truncate the Born series of the exact 
equation ([2]) to include the interaction V only in lowest 
order, i.e. 

(G) = (g) + (gVg) + ■■■ (34) 

Then one can show that, without any further approxi- 
mation, the second term of the r.h.s. can be written 

(9Vg) = (g)V(g). (35) 



Thus to lowest order in V the exact equation (j34f is iden- 
tical to the CES result Eq. (0. 

Only for intermediate coupling, where the aggregate 
H-band absorption is shifted into the region of strong 
monomer absorption, do CES and DD results not agree, 
except qualitatively under low resolution. This raises the 
question as to why the CES method applied previously 1 
produces very good agreement with measured H-band 
spectra. The answer has been shown to be that the pro- 
cedure of calculating aggregate spectra using stick spec- 
tra with only one vibrational mode pro monomer taken 
into account and then fitting to a measured continuous 
spectrum by convoluting with Gaussians, is not a good 
strategy. A good reproduction of aggregate H-band spec- 
tra is only obtained in the CES method when the ex- 
perimental monomer spectrum is used. Presumably this 
is because the DD model includes only one vibrational 
mode per monomer and does not contain any influence 
of dissipation from the environment; the eigenvalues are 
discrete. When continuous experimental monomer spec- 
tra are used, such effects are included implicitly and the 
only assumption of the CES method in its generalised 
form 19 is that such dissipative effects act in the same 
way on both monomer and aggregate spectra. 
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